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Abstract 

oo : 

The preasymptotic effects originating from the four-quark operators are disentangled from the other 
i contributions by considering appropriate combinations of inclusive decay rates. Under the assumption 
k*" . of the isospin and heavy-quark symmetry, a set of relations connecting charmed and beauty decays 
] is obtained without invoking specific models. The results are compared with other approaches and 
CN * confronted with experiment. 
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It is astonishing that a lot of physical observables, decay rates of semileptonic, nonleptonic, 
and radiative decays are described in the framework of the inverse heavy-quark mass expansion 



5_i ■ in terms of relatively few basic quantities, e.g., quark masses and hadronic expectation values 
(HEV) of several leading local operators [|I|, ||, [|. The underlying theory is based on a few 
field-theory basics, such as the operator product expansion, quark-hadron duality, and certain 
well-known symmetries. Even more surprisingly, the theory seems to work rather well in the case 



of charmed hadrons, where the expansion parameter \J^Q{D)/m 2 c ~ 0.5 is by no means very 
small. A systematic analysis leads to very clear lifetime-pattern predictions, the agreement with 
experiment being reasonable even for absolute lifetime values [[§, fj). 

In beauty hadron decays, one expects that the whole theory should work much better, since 



the expansion parameter, \j \Xq{B) j ~ 0.13 <C 1, is significantly smaller than in the case of 
charmed hadrons. Still, in spite of the overall agreement between theory and experiment [I], there 
are some questions to be answered. The main question is: does the quark-hadron duality work? 
Unfortunately, for a reliable study of duality, one needs complete control over nonperturbative 
phenomena, which is out of reach of the present theory. Nevertheless, the validity of duality was 
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studied in exactly solvable 't Hooft model in QCD 2 [§], (7|, leading to a perfect matching, which is 
certainly encouraging. Still, one would be happy to find stronger support for duality in the (3 + 1) 
theory, too. 

The question is: could one try to reduce the uncertanties of the calculation, like the strong 
sensibility to the value of ttiq, and/or find the way to disentangle various preasymptotic effects 
which cause lifetime differences? 

The decay rate of a decaying hadron is generally of the type 
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where c{) are the Wilson coefficients and (Hq\Od\Hq) are the matrix elements of the D-dimensional 
operators which are suppressed by the inverse power of mass I/wq -3 . 
The following comments are in order: 

1. It would be desirable to extract a combination of rates T' that depends only on the operators 
with D = 6, since they are the operators responsible for lifetime differences. 

2. Since the D = 6 operators are suppressed by m^ 3 with respect to the leading D = 3 operator, 
this would reduce the mass dependence of T' to ttiq, thus significantly reducing the errors 
coming from the uncertainty of vtiq. 

3. Even truncating the sum in (]T|) at D = 6, one would presumably be able to test the role of 
preasymptotic effects with satisfactory accuracy, relying on heavy-quark symmetry. That is 
the aim of this paper. 

As both the decays in the b and c sectors of hadrons are described by the same formalism, 
it seems worthwhile to investigate the possibility of establishing connections between these two 
sectors. In a recent paper M, Voloshin demonstrated one of such possible connections. The 
strength of that approach is in avoiding the model dependence assuming SU (3) flavor and heavy- 
quark symmetry (HQS). One can extend this approach to obtain more predictivity. Using the 
matrix elements extracted from charmed baryons, a parameter Fg (parametrizing four-quark 
contributions to beauty baryons) can be determined and concrete numerical predictions can be 
obtained in this way ||. This procedure brings the ratio r(A(,)/r(S^) to better agreement with 
experiment and gives predictions for the lifetimes of beauty baryons. 

In the present paper we adopt a similar, yet different strategy from that followed in ||. In ||, 
the author, when establishing the connection between c and b baryons, always relates baryons with 
the same light-quark content, i.e., mutually related baryons differ in the flavor of the heavy quark. 
Also, the author explicitly extracts values of four-quark operator matrix elements and then applies 
them elsewhere. Our approach differs in both these respects. The heavy hadrons that we relate 
neither contain the same heavy quark(s) nor the same light (anti) quarks - they are affected by the 
same type of the four-quark operator contribution. Furthermore, we form such combinations that 
the four-quark operator matrix elements get reduced. 

In forming such combinations, we consider the leading modes of the decay of c and b quarks 
with respect to the CKM matrix elements. For the c quark, we then have only one nonleptonic 
mode, c — > sdu, and one semileptonic mode per lepton family, c — > slvi. In the case of b quark, 
there are two nonleptonic modes, b — > cud and b — > ccs, as well as one semileptonic mode per 
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lepton family, b —>■ clv[. As we consider b baryons containing only light quarks along with b 
quarks, the semileptonic Cabibbo leading modes for the decay of the b quark do not appear. We 
also disregard mass corrections of four-quark operators coming from the massive particles in final 
decay states. This is a better approximation for c decays (as m 2 s /m 2 c ~ 0.01) than for b decays 
(where m 2 /m 2 ~ 0.1). 

Traditionally, the effects of four-quark operators are called the positive Pauli interference, the 



negative Pauli interference and, W exchange or annihilation |nj, [TTJ. We consider pairs of one 
charmed and one beauty particle which contain the same type of four-quark contribution. It is 
important to notice that in b and c decays different light-quark flavors participate in the same type 
of four-quark contributions. Using isospin (SU(2) flavor) symmetry and HQS we are able to form 
appropriate decay-rate differences which lead us to the final result. 

We start our considerations with heavy mesons. The first pair of considered particles are 
D + (cd) and B~ (bu). In both these mesons the effect of negative Pauli interference occurs. The 
contributions of four-quark operators to the decay rates of these two particles are described by the 
same type of operators. Moreover, the operators for the b case can be obtained from those in the 
c case by making the substitution c — > b, d —> u. The second pair of particles form D° (cu) and 
B° (bd). In both these mesons the effect of W exchange occurs. 

Next, we would like to isolate the effects of four-quark operators, i.e., we have to eliminate 
the contributions of the operators of dimensions 5 and lower. We achieve this goal by taking 
differences of the decay rates of D and B mesons assuming the isospin symmetry. This procedure 
leads us to the following relations: 

T(D+) - T(D°) = ^A\ Vcs \ 2 \V ud \ 2 [(D + \P cd \D + ) - (D°\P CU \D )] , (2) 
n 2 m 2 

T(B-) - T(B°) = -J^\V cb \ 2 \V ud \ 2 [(B-\P bu \B-) - (B°\P bd \B )] , (3) 

where P's denote the appropriate four-quark terms from the heavy-quark effective Lagrangian, 
see, e.g., Ref. ||. 

Using HQS and isospin symmetry we can express the terms describing the negative Pauli 
interference as 



P cd = E <?0*)O?(n) = E ctl 9 MO% g M + 0(l/m c ) , (4) 

i=l i=l 



P bu = E PMopQ*) = E ^,L 9 Mo$ eg M + o(i/m b ) , (5) 

i=l i=l 

where Q denotes a heavy quark in the heavy-quark limit mg — > oo. For the matrix elements 
of the above mentioned operators we obtain (by applying HQS and isospin symmetry again) 
(D+\P cd \D+) = A negint + 0'(l/m c ) and (B~\P bu \B~) = A negint + 0'(l/m b ), where the matrix 

element is A negint = (Mq 5 | Ya=\ c?,negint{^)0®negint{^)\ M Qq) an d M Qq denotes the mesonic state 
with one heavy quark Q and one light antiquark q in the heavy-quark limit. In the preceding 
relations we have also stated that heavy-quark mass suppressed corrections to the operators and 
their matrix elements do not have to be the same. 
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In an analogous manner we can obtain the expressions for the operators describing W exchange 
using the notation explained above: 



p cu = E <»ofV) = E $LM°?,LM + 0{ M m <) > ( 6 ) 

i=l i=l 

p m = E 4 d (M d (ri = E ^Mo?LM + o(i/m b ) . (7) 

8=1 i=l 

The matrix elements of the operators given by (§) and (0) are (D°|P CU |D°) = A exc h + 0'(l/m c ) 
and (fi°|P M |£°> = A exch + 0'(l/m b ), with A exc/l = {M m \ E* =1 cg*^) OgLhMI M Q5 ). In all 
relations given above, the dependence of the coefficients q and the operators Oj on the scale \x is 
explicitly displayed. Since the dependences of the operators and the coefficients on /i cancel each 
other, there is no need for an explicit value for the /i scale, and also no hybrid renomalization is 
required. That scale, however, should be small compared with the heavy-quark masses, to allow 
for the proper heavy-quark expansion and the reduction of the four-quark operator contribution 
in the heavy-quark limit, as shown below. 

Combining all preceding expressions, we obtain for the decay rate differences 

G 2 m 2 

Y{D + ) - T(D°) = -J-^\V cs \ 2 \V ud \ 2 [A negmt - A exch + 0(l/m c )} , (8) 

r(£T) - T(B°) = 9jfi\ Vcb \i\ Vud \2 [A _ Aexch + 0(l/m b )} . (9) 

As m b 3> m c , it is clear that the approximation made above should work much better for beauty 
particles than for charmed particles. Still, we assume that the approximations made are justified 
in both c and b decays and that the corrections of order 0{l/m c ^) are not large. In this case, we 
obtain the final relation for mesons: 

Rn T(B-)-T(B°) m 2 \V cb \ 2 

r s rM - iV) = + ° (llm " 1/mt) ■ (10) 



Let us consider relation ( |10|) in more detail. By forming the lifetime differences (0) and 
(§) we have not only eliminated the effects of the leading operators in the mesonic decay rates, 
but applied the HQS first at the subleading level of 0(l/m z cb ). This has in turn enabled us to 



eliminate the dependence on four-quark matrix elements in relation fllQ|) . Also, the sensitivity to 
the choice of heavy quark masses is now significantly reduced since in expression (PH|) we have 
only the second power of masses. The procedure, however, has its limits. It enables us to reduce 
the matrix elements, but cannot be extended to CKM suppressed modes. Also, mass corrections 
due to massive particles in the final decay states spoil the procedure and, therefore, have been 
left out. Nevertheless, both of these corrections are under good theoretical control. Moreover, the 
contributions from mass corrections and suppressed modes are smaller by more than one order of 
magnitude than the leading contribution and cannot significantly change relation (fLOf). 

Starting from expression (0), we may check the standard formalism of inclusive decays against 
experimental data, especially if the four-quark operator contributions are sufficient to explain the 
lifetime differences of heavy mesons. To perform such a check, we determine the quantity r BD using 
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the experimental values for lifetimes and theoretically (using relation flit])) and then compare the 
values. Taking the experimental values from ||12|| , we obtain 

= 0.030 ±0.011. (11) 

The large error of rfj^ comes from the fact that the experimental errors of lifetime measurements 
of B~ and B° are large compared with the difference between central values of the results for B~ 
and B°. Inspection of |12j shows that the experimental situation in the sector of B mesons is far 
from being settled and more precise measurements of B meson lifetimes are needed to reduce the 
uncertainty in the value of rfj^. 

The theoretical value rf h D is calculated from expression (|10|) . The numerical values for heavy- 
quark masses are taken to be m c (m c ) = 1.25 ± 0.1 GeV [[HJ and rrib(lGeV) = 4.59 ± 0.08 GeV 
|i"4fl for the reasons discussed in || . The values for the matrix elements of the CKM matrix are 



taken from |T_2[ to be |V CS | = 1.04 ±0.16 and \ V c b\ = 0.0402 ± 0.0019. Using these numerical values 
we obtain 

rf h D = 0.020 ±0.007. (12) 



The relatively large error in (12) comes predominantly from the large experimental error of the 
| V cs | CKM matrix element . 

Direct comparison of the numerical results ([□]) and (|T2"D indicates that these are consistent 
within errors. Since the relative errors of both results are rather large and experimental data are 
still fluid, it is possible that these results will experience some change. Nevertheless, we expect 
that the two results will remain comparable and even closer to each other. 

These results confirm, in a model-independent way, that four-quark operators can account for 
the greatest part of decay rate differences of heavy mesons, i.e., that the contributions of other 
operators of dimension 6 and higher-dimensional operators are not so important. Of course, these 
conclusions should be interpreted in the light of the approximations made. 

The procedure presented for heavy mesons can be extended to heavy baryons. First, we 
consider the system of singly heavy baryons. These baryons contain two light quarks, which 
introduces two types of four-quark operator contributions for each baryon. Again, we choose 
two pairs of singly heavy baryons, each pair containing one charmed and one beauty baryon 
related by the same dominant type of four-quark contributions. The first pair contains (cus) 
and (bds), while the second pair contains S° (cds) and S° (bus). The first pair exhibits the 
negative Pauli interference and different nonleptonic and semileptonic contributions of four-quark 
operators containing s quark fields. The second pair comprises the 1^-exchange effect and the same 
contributions of four-quark operators containing s quark fields as in the first pair of baryons. It is 
important to notice that the four-quark operators that contain s quark fields do not describe the 
same type of effects in c and b decays in our pairs (and therefore are not given by the contributions 
of the same type in the effective Lagrangian). Nevertheless, we form decay-rate differences in 
such a manner that the contributions of four-quark operators containing the s-quark field cancel 
(assuming isospin symmetry). This cancellation is important in the case of b decays because in 
the contributions of s-quark type four-quark operators significant mass corrections appear because 
of two massive c quarks in the final state. 

In the singly heavy sector we, therefore, form the following differences: 

r(s c + ) - r(H°) = ^|^| 2 |K,| 2 [(H+|P c 1~+> - (H°|P C 15°)] , (13) 



5 



r(sr) - r(sg) = ^l^l 2 IK,l 2 [(s 6 -|P M |s,-) - <H°|P b ls°>] . (w) 

Following the procedure displayed in relations ([|) to fllCf ) and using the isospin symmetry 
and HQS again, we obtain the expressions (note that the four-quark operators which contribute 
to the negative Pauli interference in mesons contribute to W exchange in baryons and vice versa 



mil) 



and 



r(H+) - r(S°) = 9jfZ\ Vc f\V ud \ 2 [B negint - B exch + 0{l/m c )\ , (15) 
r(Sr) - r(Sg) = ^^\V cb \ 2 \V ud \ 2 [B negmt - B exch + 0{l/m h )\ , (16) 



An 

where the matrix elements of the four-quark operators between singly heavy baryons Bq qq / have the 

iormB negint = (B Qqql \ T^^^neginti^Oi^eginti^l^Qqq') arici B exch = (B Qqq >\ Ei=l C i,Lch(^)°i,Lch(^) I B Qqq') ■ 

The final relation for singly heavy baryons looks like 

r ^W^S) = ^ li+o{iMmi)l (17) 

This relation can be used to test the present model-dependent calculation of heavy-baryon 
lifetimes. Lifetime splittings both in the charmed and in the beauty sector have been investigated 
in several papers 0, 0, [|, |9], [15]]. We have recalculated our predictions from |4| for charmed 
baryon lifetimes and from || for beauty baryon lifetimes with the numerical parameters preferred 
in this paper and have used the recalculated predictions to test the relation (0), which we express 

in units ^ L 2 thus reducing the dependence on heavy-quark masses. In the case when the same 
approximations (neglecting mass corrections and Cabibbo-suppressed modes) are made, the above 
relation differs from unity by —12%. This number shows the deviation of the model-dependent 
calculation (where the contributions of four-quark operators are explicitly evaluated) from the 
model-independent prediction given only by the ratios of heavy-quark masses and Cabibbo matrix 
elements. The complete calculation with the mass corrections and Cabibbo-suppressed modes 
gives 0.79 for the above ratio, which indicates the order of neglected corrections to be less than 
10%. 

There are several important implications of relation ([H])- We can obtain a model-independent 
prediction for the still unmeasured lifetime difference between singly beauty baryons in the isospin 
doublet. Using the same numerical values for the parameters entering the right-hand side as before, 
and taking the experimentally measured lifetimes of singly-charmed baryons and S° from [I2|, 
we obtain 

r(H 6 -)-r(H°) = -(0.14 ±0.06) ps" 1 . (18) 

This lifetime difference can be compared with the model-independent prediction obtained 
by using the HQS and SU(3) symmetry in |§, T(S^) — T(S°) = —(0.11 ± 0.03)ps _1 , and some 
moderate model-dependent prediction from ||, T(S^) — r(S°) = — 0.094 ps -1 . It is also worth 
mentioning that owing to the neglected mass corrections in final decay states of beauty baryons in 
the derivation of the ratios ( |TTD and ([18]), there is no splitting between the nonleptonic rates of H° 
and Ab- Since the Cabibbo-suppressed modes have also been neglected , the total decay rates of 
S° and Ab appear to be equal at this level. Therefore, the predictions from (|T7]) are also valid for 
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the r(H£") — r(Aft) difference. We can see that the results for the predicted splitting in the isospin 
doublet of beauty baryons obtained using different methods are all consistent. 

It is interesting to examine the explicit expression taken from ||, eq.(23), which can be 
rewritten as 



r(sr) - r(s°) = 4^ [°- 9ir ( s c + ) - o.85r(s°) - o.oer(A c )] , (19) 

mi Vr 1 



and which is obtained by applying £77 (3) symmetry, and with Cabibbo subleading effects included 
in charmed baryon decays. The appearance of the total decay rate for A c is the direct consequence 
of the SU(3) relations used. Otherwise, we can note a very similar structure to our relation 
(|17l) obtained within the £77(2) approximation. The first two coefficients in front of the rates in 
(PH) differ from one by less than 15%, where about 5 — 10% difference comes from the Cabibbo- 
suppressed modes included, and the rest comes presumably from the difference between the £77(3) 
and SU (2) approximations applied. 

In the last case, we apply our procedure to the system of doubly heavy baryons. Again, we form 
pairs of one doubly charmed and one doubly beauty baryon. The first pair includes 5+ + (ecu) and 
S bb (bbd). The decays of both baryons in this pair include the effect of negative Pauli interference. 
The second pair is formed from (ccd) and S bb (bbu). Both of these baryons exhibit the effect 
of W exchange. There are no Cabibbo-leading semileptonic four-quark contributions in doubly 
heavy baryons. Appropriate decay-rate differences are given by the following expressions: 

r(H++) - r(~+) = ^|K s | 2 |K,l 2 [(Hi + |P c lH+ + ) - <s+|Pi ~+ c >] , (20) 
r(3£) - n~° bb ) = ^\v^\vM{^ b \P u \^ b ) - £ bb \P bu \z bb )] ■ (21) 

Using isospin symmetry and HQS, we obtain, per analogiam with above derivations, relations 
for the decay-rate differences of doubly heavy baryons which lead to the final relation for doubly 
heavy baryons 

rS) - r(=gj) ml\v cl ' 2 



■ -r (S . + )-r( 3a ^[i + oa/^i/™,)]. W 

Again, we have achieved the reduction of the matrix elements of four-quark operators between 
doubly heavy baryons in the heavy quark limit. 

Relation (|22| ) enables us to estimate the difference in the decay rates of doubly-beauty baryons 
using recent results for doubly-charmed baryon lifetimes. Using the calculated values T(S+ + ) = 
0.952 ps^ 1 and T(S+) = bps" 1 with the values of the parameters m c = 1.35 GeV, m b = 4.7 GeV, 
\V cb \ = 0.04 and \V CS \ = 1.04 gives 

T(~ bb )-T(~° bb ) = -0.073 ps- 1 . (23) 

So far we have applied our reduction procedure separately to heavy mesons, singly-heavy 
baryons, and doubly- heavy baryons. Still, by inspection of relations (|T0|), flTT|), and (|22| ) we see 
that r BD , r bc , and r bbcc are the same up to 0(l/m c , l/rrib), i.e., 



r(B-) - r(s°) r(s 6 -) - r(sg) r(E bb ) - r(s° 6 ) ml \v c 



d>\ 
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r(D+) - r(D°) r(s+)-r(sg) r(H++)-r(s+) m *\v c 



(24) 
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Thus, we obtain a relation in (|24] ) which clearly indicates certain universal behavior in the de- 
cays of all heavy hadrons. The existence of some universality could have been anticipated from the 
fact that the same expression (|XJ) describes the decays of all heavy baryons. Still, that universality 
attains its concrete, model- independent form in the relation fl2"4"|) . This relation connects all sectors 
of heavy hadrons that are usually treated separately: mesonic and baryonic, charmed and beauty. 
Also, this relation brings some order in the otherwise rather intricate pattern of heavy-hadron 
lifetimes. An advantage of this relation is that by knowing some decay rates, one can calculate or 
give constraints on some other decay rates. Also, knowing decay rates from experiment, one can 
test the findings of the theory in a model-independent fashion. These results will help to establish 
the limitations of the present standard method of calculating inclusive processes and to test its 
underlying assumptions. 
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